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Abstract 



We study the physical electron in quantum electrodynamics expanded on 
the light-cone Fock space in order to address two problems: (1) the physics of 
the electron's anomalous magnetic moment Og in nonperturbative QED, and 
(2) the practical problems of ultraviolet regularization and renormalization 
in truncated nonperturbative light-cone Hamiltonian theory. We present re- 
sults for Oe computed in a light-cone gauge Fock space truncated to include 
one bare electron and at most two photons; i.e., up to two photons in flight. 
The calculational scheme uses an invariant mass cutoff, discretized light-cone 
quantization (DLCQ), a Tamm-Dancoff truncation of the Fock space, and a 
photon mass regulator. We introduce new weighting methods which greatly 
improve convergence to the continuum within DLCQ. Nonperturbative renor- 
malization of the coupling and electron mass are carried out, and a limit on 
the magnitude of the effective physical coupling strength is computed. A 
large renormalized coupling strength an = 0.1 is then used to make the non- 
perturbative effects in the electron anomalous moment from the one-electron, 
two-photon Fock state sector numerically detectable. 
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I. INTRODUCTION 



Many years ago Feynman issued the following challenge [|T|: "It seems that very little 
physical intuition has yet been developed in this subject [of quantum electrodynamics]. In 
nearly every case we are reduced to computing exactly the coefficient of some specific term. 
We have no way to get a general idea of the result to be expected. ... As a specific 
challenge, is there any method of computing the anomalous moment of the electron which, 
on first rough approximation, gives a rough approximation to the a term and a crude 
one to o?] and when improved, increases the accuracy of the o? term, yielding a rough 
estimate to o? and beyond?" This challenge was taken up by Drell and Pagels 0, who 
used a sideways dispersion relation and low-energy theorems for Compton scattering [|| to 
construct consistency conditions for the anomalous moment. Their approach did meet with 
some success, particularly in understanding the sign of the a term; however, the dispersion 
relation requires an ultraviolet cutoff, and low-energy approximations of the integrand are 
not completely adequate. 

We propose to meet Feynman's challenge by using discretized light-cone quantization 
(DLCQ) 0. By constructing the dressed electron state in Fock space we can in principle 
compute physical properties of the electron exactly In practice, various truncations 

are required, but the approach remains nonperturbative. Instead of producing an expansion 
in a, we produce an expansion in Fock particles, i.e. the number of photons in fiight. The 
computation is equivalent to a selective summation of graphs to all orders, but is actu- 
ally done by diagonalizing a matrix approximation of the light-cone mass-squared operator. 
In this form the calculation becomes a testing ground for techniques of nonperturbative 
renormalization. 

When two-photon intermediate states are allowed, graphs such as the multiloop graph 
in Fig. |l|(a) enter the calculation and are summed to all orders. Even at the one-photon 
level, the calculation is nonperturbative because there are infinite-order contributions from 
graphs of the sort in Fig. |l](b) and (c). However, in the case of Fig. |l](c), crossed-photon 
graphs and Z graphs are needed to cancel a divergence at zero longitudinal momentum for 
any instantaneous electron. Because of these cancellations, we place the nonperturbative 
part of the one-photon contribution into a two-photon calculation. 

The selection of the graphs to be summed is driven by the truncations made and is 
designed to make the calculation tractable. The truncation in particle number is physically 
reasonable; the work of Drell and Pagels |0] shows that states with few photons in fiight 
are dominant. As the actual number of particles is varied and as different truncations of 
the interactions are explored, one can gain a better understanding of the physics of the 
anomalous moment. Thus our calculation can be viewed as the beginning of a possible 
program, with systematic improvements available. It might even become a model for how 
to proceed with nonperturbative calculations in quantum chromodynamics (QCD). 

In the work presented here we regulate the theory by an invariant-mass cutoff, which 
limits the light-cone energy of the Fock states included in the basis, and by a Tamm- 
Dancoff truncation |^ of the number of constituents. These restrictions keep the numerical 
calculation to a very reasonable size but they do complicate the renormalization. The 
truncations turn the bare parameters of the theory into Fock-sector dependent functions of 
momentum and require careful construction of appropriate counterterms [|ll| to correctly 
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approximate the solution to the original theory [0. We fix these functions by applying 



conditions on mass eigenvalues and vertices in the presence of spectator constituents. 

Two alternative renormalization procedures now exist. One is the similarity transforma- 
tion developed by Wilson and Glazek and Wegner |jT^, where counterterms are generated 
perturbatively as the Hamiltonian matrix is narrowed in the range of allowed light-cone en- 
ergies; the final Hamiltonian matrix is used nonperturbatively. This approach has been 
applied by Perry and co-workers [|T^]. The other is based on introduction of Pauli-Villars 
regulators [|16| before the quantization and numerical schemes are selected, so that coun- 
terterms can be supplied simply by adjusting the bare parameters. This approach has just 



recently been successfully tested in a simple model |T^, and it should soon be considered 
for the anomalous-moment problem studied here. 

We compute the anomalous moment [^8|,|19| from a spin-fiip matrix element of the plus 



component of the electromagnetic current pO[- We approximate the Fock-state expansion of 



the dressed electron with a truncation to no more than two photons and one electron. The 
eigenvalue problem for the wave functions becomes a coupled set of three integral equations. 



To construct these equations we use the light-cone Hamiltonian derived by Tang et al. [21 



regulated by the invariant-mass cutoff. The photon mass is taken to be one tenth of the 
electron mass, to help control infrared divergences. The coupling strength is set at a = 1/10, 
because limitations on numerical accuracy make nonperturbative effects discernible only at 
large coupling. The calculation is not an attempt to compete with the accuracy of the 
perturbative calculations by Kinoshita and co-workers p^ . 

The bare electron mass in the one-photon sector is computed from the one-loop correction 
allowed by the two-photon states, where one photon is a spectator. We then require that 
the bare mass in the no-photon sector be such that physical mass is an eigenvalue of the 
light-cone Hamiltonian. 

The three-point bare coupling is also sector dependent. There are no vacuum polarization 
effects, because pair production is removed by the Tamm-Dancoff truncation. However, the 
truncation violates the Ward identity so that vertex and wave function renormalization do 



not cancel |2^. A consequence of this is that the physical coupling is limited to a cutoff- 
dependent finite range of allowed values. We compute the critical coupling that defines the 
upper limit. 

The vertex renormalization is fixed by considering the proper part of the transition am- 
plitude for photon absorption by an electron at zero photon momentum. A means by which 
this transition amplitude can be computed from the lowest eigenstate is constructed. Full 
diagonalization of the Hamiltonian is not required; however, the renormalization condition 
and the eigenvalue problem must be solved as a coupled system. The wave function renor- 
malization is directly available from the bare amplitude in the Fock state expansion. 

To within the accuracy of the calculation, the values for the anomalous moment become 
constant for an ultraviolet cutoff sufficiently large. However, most four-point graphs that 
arise in the bound-state problem are log divergent. To any order the divergences cancel 
if all graphs are included, but the Tamm-Dancoff truncation spoils this. The resulting 
logarithmic effects are not detectable in the numerical results. 

The calculations presented here build on the significant amount of work that has been 
done recently on the use of light-cone quantization [^,^ p4tp5| in the construction of solv- 
able bound-state problems for strongly interacting theories. The coordinates used are based 
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on the choice of t + z as the time coordinate, where t is the ordinary time and z any 
Cartesian spatial coordinate. A variety of field theories have been analyzed in this way 
J6|, pT] , p7| -|41[] . Most theories considered have been simple model theories in one space 



dimension |]6|j26|- |33| , ^^ ; however, there have been studies of three dimensional theories, 
including the Wick-Cutkosky model P2| , p^ , |35| , the Yukawa model [^], quantum electro- 
dynamics (QED) [^p7| , and quantum chromodynamics (QCD) [^. There has also been 
some work on nonperturbative scattering calculations P5|-^] and on the stationary phase 
approximation to the soliton in 



Much of this work has involved numerical studies. Brodsky, Pauli, and co-workers have 
analyzed various one-dimensional theories (Yukawa 0, QED ||2^, and QCD |^) and have 
devoted a considerable amount of effort to three-dimensional QED |pll , |37| . Some work on 
the three-dimensional Wick-Cutkosky model |^ has been done by Sawicki and co-workers, 
Ji and Furnstahl and Wivoda and Hiller One-dimensional scalar theories, 0^ 

and 0^, have been studied by Harindranath and Vary Work on the one-dimensional 
Yukawa model has been done by Harindranath, Shigemitsu and Perry this was based on 
a Tamm-Dancoff truncation and used basis-function methods as well as a discretization 
technique. The basis-function methods have been extended to the three-dimensional case by 
Glazek et al. A preliminary treatment of QCD in three dimensions has been attempted 
by Hollenberg and co-workers [Q. Dimensional reduction of QCD to an effective theory in 
1 + 1 dimensions has also been considered ||40|j4l[| . 

Other nonperturbative approaches applicable to QCD include lattice gauge theory 
47| , |48[| , sum rules |^9|, and Schwinger-Dyson equations DSOf. A particularly successful form 



of lattice theory has been developed by Lepage and collaborators |^ who use tadpole- 
improved actions to reduce discretization errors. The Hamiltonian form of lattice theory 
I p^ is actually similar to the approach usually taken in light-cone quantization, in that a 
Hamiltonian operator is constructed and partially diagonalized [Q. There has also been 
work on combinations of the lattice with light-cone methods in the transverse lattice method 
1 5^ and in direct use of a light-cone lattice [^ . 

Two important aspects of QCD that all these methods address are vacuum structure 
and symmetry breaking. In light-cone quantization, the vacuum appears to be the trivial 
perturbative vacuum. This has the advantage that one can compute massive states imme- 
diately without first computing the vacuum state. However, in equal-time quantization, the 
nontrivial structure of the QCD vacuum is known to be important. This paradox of the 
trivial vacuum has received much attention. Nonperturbative analyses of various light-cone 
models indicate that interactions induced by zero modes ||56|- p^] and other considerations 
]TT| play important roles in generating effects such as symmetry breaking |^,^| that are 
usually associated with the vacuum. 

The progress made recently in the application of light-cone quantization owes much to 
earlier work ^1| . The development then was aimed at perturbation theory, and in particular 
its application to deep inelastic scattering, but much has been carried over to bound-state 

m. 



problems. New work on perturbation theory has also been done 

An outline of the remainder of the paper is as follows. The discretized light-cone formu- 
lation of the anomalous moment problem is given in Sec. 0. The nonperturbative mass and 
coupling renormalization that we use are described in Sec. |T|. Finite corrections associated 
with photon zero modes and with ambiguities in infinite renormalizations are discussed in 
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Sec. Our results are presented in Sec. 0, and a brief summary is given in Sec. ^ 



II. DISCRETIZED LIGHT-CONE QUANTIZATION 
A. Light-cone quantization 

We define light-cone coordinates by 

x^=t±z, x± = {x,y). (2.1) 

Momentum variables are similarly constructed as 

p^ = E±p,, P± = {Px,Py)- (2.2) 

The time variable is taken to be x~^, and time evolution of a system is then determined by 
V~ , the operator associated with the momentum component conjugate to x~^. Usually one 
seeks stationary states obtained as eigenstates of V~ . Frequently the eigenvalue problem is 
expressed in terms of a light-cone Hamiltonian (mass-squared operator) 

i^Lc = V+V- - rl (2.3) 



as 



i^LC^ = M^^ , (2.4) 

where M is the mass of the state, and and 'P± are momentum operators conjugate to 
x~ and x±. 

It is convenient to work in a Fock basis {\n : pt,p±i)} where and 'Px are diagonal, 
with n the number of particles and i ranging between 1 and n. To simplify the notation 
only one particle type is included explicitly. The state \1/ is given by an expansion 

^ = E /H" [d^k^UM^^ k±)\n : p+ = xP+,px = xP^ + k^) , (2.5) 

with 

n n J n n >2l . 

[dx]. = 4vr5(l-E^.)n/%^ [rf^fcJ„ = 4vr^5(Efc^,)nY#, ^^.G) 

1=1 i=i ^^v^« i=i i=i 

(P^, P_l) the total light-cone momentum, and ipn interpreted as the wave function of the 
contribution from states with n particles. The solution of (p.4|) in principle yields these wave 
functions. 

The anomalous moment a^. is computed from the standard form factor -^2(9^) at zero 
momentum transfer: 

ae = i^2(0) . (2.7) 



In the standard light-cone frame ||63] where 
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g = (0,gi/p+,q± = gix), (2.8) 

the form factor can be computed from the spin-flip matrix element of the plus component 
of the current: 

?1 r. / 2x 1 



F2iq') = ^{p + q,UJ^mPA). (2.9) 



2me 2p 



Brodsky and Drell [20] have given a reduction of this matrix element to a convenient form 



that depends directly on the wave functions. From this we have 

[dx]n [(fk±\ni'nt{x, k±) ^ Xi——i)nlix , k±) , (2-10) 

where Cj is the fractional charge of the struck constituent. 

Up to this point, we have used formally exact expressions. A key approximation to be 
made is the truncation of all sums to a finite number of particles. The result is the light- 
cone equivalent of the Tamm-Dancoff approximation [Q. The eigenvalue problem becomes 
a finite set of equations that are in principle solvable. However, the truncation has many 
consequences for the renormalization of the theory [|T0| and for comparisons to Feynman 



perturbation theory Some of these consequences are discussed in Sec. |n. 

In addition, QED requires regularization and renormalization. To regularize it, we use a 
cutoff on the invariant mass of the allowed Fock states 

^m[+fcL<^2_ (2.11) 

This limits the relative transverse momentum k±_ of each constituent and keeps the longi- 
tudinal momentum away from zero. The latter aspect is important for control of spurious 
infrared singularities, which are discussed in |III C|. An additional cutoff, that limits the 



change in invariant mass across any matrix element of the Hamiltonian could be con- 
sidered. 

When only states with at most one photon and no pairs are retained, and instantaneous 
interactions are neglected, Brodsky and Drell |^ have shown that Eq. ( p. 101 ) reduces to 



ami r g(A^ - (m^ + A:i)/(1 - x) - (m^ + A;i)/x) 

a, J axa _ (^2 + p^y^i _ 3,) _ (^2 + p^y^^2 ' y^-^^) 

which in the limit of A — > 00 becomes 

«e = af = — / , , , • (2.13) 



2t[ Jo x"^ + (1 — x){m.y/me) 



Because the instantaneous interactions are higher order in a, this is the leading perturbative 
result. The integrals involved can all be done analytically even for finite cutoff, although the 
final form is not instructive. For = 0, the resulting formula yields the standard Schwinger 
T8| contribution of a/27r at infinite cutoff. In general this provides a point of comparison for 



numerical calculations with one or more photons. The inclusion of the dependence on the 
photon mass in the analytic result is crucial for comparison with numerical results calculated 
with nonzero m-y because the mass dependence is quite strong, as can be seen in Fig. 0. 
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B. Discretization 



The most systematic approach to discretization of the eigenvalue problem is the method 
originally suggested by Pauli and Brodsky 0, discretized light-cone quantization (DLCQ). 
In essence it is the replacement of integrals by trapezoidal approximations, with equally- 
spaced intervals in the longitudinal and transverse momenta 

p-^^-n, [—n^,—ny) . (2.14) 

The length scales L and determine the resolution of the calculation. Because the plus 
component of momentum is always positive, the limit L —>■ oo can be exchanged for a limit 
in terms of the integer resolution ^ 

K = -P+. (2.15) 

TT 

The combination of momentum components that defines i^LC is then independent of L. The 
longitudinal momentum fractions Xi become ratios of integers rii/K. Because the Ui are 
all positive, DLCQ automatically limits the number of particles to no more than K. The 
integers Ux and ny range between limits associated with some maximum integer A^_l fixed 
by the invariant-mass cutoff. A finite matrix problem is then obtained without an explicit 
Tamm-Dancoff truncation; however, this number of particles is much too large in practice 
for numerical treatments of three-dimensional theories. 

We use antiperiodic boundary conditions for the fermions and periodic boundary condi- 
tions for the photons. These restrict the integers n associated with longitudinal momenta to 
being odd for fermions and even for photons. The description of the dressed electron state 
must then use odd values of K. 

In most applications, DLCQ is introduced at the level of second quantization. This can 
yield a compact expression of the eigenvalue problem. Recently, a transformation of the 
DLCQ Hamiltonian to a Gaussian basis has been suggested [^]; however, the steps for 
renormalization in that basis have not been worked out. 



The application of DLCQ to QED is summarized in pT|, which we use as a starting 



point. This includes use of light-cone gauge [Q, with = 0. Modifications of this gauge 



choice due to zero modes ||5j are discussed in Sec. 

The fundamental interactions of light-cone QED are illustrated in Fig. |^. For the cal- 
culation reported here we do not include any pair production processes. The instantaneous 
photon interactions are then completely excluded because each Fock state has only one 
fermion. We also exclude the fourth diagram (and its conjugate) to decouple two-photon 
states from the bare electron state; this simplifies the calculation and limits the role of 
two-photon states to that of providing the basis for inclusion of crossed-photon graphs. 
There is also a technical modification of the interaction associated with the third diagram of 
Fig. ^ which is discussed in Sec. [Ill C| . After the exclusions have been made, the light-cone 
Hamiltonian becomes 

ml + {n_Ln/L±f . , ml + {m_i7i/L_Lf . 

HlC = T. E '/r^ -blsbn,s + T. E ZJl? <,^<^rn,^ (2-16) 

n s=±l/2 'y ^ m A=±l '"'/ ^ 
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FIG. 3. Complete set of diagrams for the fundamental interactions of QED in light-cone 
quantization. Solid lines represent fermions; wavy lines represent photons. A line with a bar 
through it indicates an instantaneous interaction expressed in the Hamiltonian as a four-body 
operator. Only the first three diagrams are included in the present calculation, and for the third, 
the piece kinematically equivalent to a Z graph is neglected. 
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I- 



m n2 



.=±^/2 V^l"^2 Ui + mi -2' -2' 



1112 — ni 



with e^^A = — (A,i)/\/2. 

The discrete form of the spin-s eigenstate is 



I»,Zn si,Ai 

+ ^n+rn,+rn^,K Y V'2.(zi, ZHi, 1^2; «!, Ai, A2) ^6^,,,a^^^;,^a^^^;,j0) , (2.17) 

where X = [K, K± — 0). According to the eigenvalue equation i?LC^ = M^^ the amph- 
tudes ipis must satisfy the following (discretized) integral equations: 

(M^ - mD^os = eo ^^^r (^n+m.i^ ^= f - - ^) ^is{n, m; -s, 2s) (2.18) 



+ e± _2s -2s) > , 

m ) 



ml{n') + (n'_L7r/L_L)2 _^ "^7 + {m'^'^/Lj 



n'/K m'/K 

K 



^is(n',m';si, Ai) 



47rL 



-ei(2i')ei(2l)(^Ai,2si^^^^^ — ^V'is(21, m Si, 2si) 
n — m' 



Knif. 1 /I 1 \ 
+^o^rW — n \ ^\i,2s^si-s^os (2.19) 

/Ta/tt if* / rn'j_ n'j_\ ^ m'^ c \ , 



2V7r-L_L^ y/2m\n n + 



m 
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X {^2s{n,m,m'; -Si,2si, Ai) + il>2s{n, mf , m; -Si, Ai,2si)} 



1 [ 

V2L1 r;^. V2m I 



V^-^i n,_ 



2si 



m n 



m n 



'ip2s{n,r]l,m-, si, \i,2si) 





_< 


m 


n' 


m± 




m 


n' 



;- 4^2s{n,m^,m; Si,Xi, -2si) 



and 



n'/K 



+ 



m'jK 
ml + (m'_L2vr/L_ 



\2-\ 



m'2/K 



i'2s{n,m!.i:ml2^ si, Ai, A2) 



KrUp 



2V27rL 



JA2 ,-2sr 



ei(?i' + rn/2) ( 1 



m'2 n' + m^ 



+^Ai,-2si 



61(22' + m'l) / 1 



1 



Ky/ri i ei(n' + 772,2 



(5. 

+'^Ai,2si 
+5a2 -2si 



.2; ,* 



ra' + nil 



m^2 ^ 



7 + mi'i,?2i2; -■si^ -^2) 



m2 

61(22' + 22l'i) 



A2 ■ — 7 V^i.(li' + 22l2,22l'i;si,Ai) 



7712 ^ 

'm'_i_i n'_i_ 



(2.20) 



^±,Ai 



mi 

61(22' + 22I2) * 



7 T V^is(l^' + 22l'i'12l2; -515-^2) 



mi 72 
^±2 



m' o n'l + m' 



mb 



ra' + m 



— ^is(li' + 22l2,I?l'i;si,Ai' 



-2si " 



61(22' + Wll) * / ''^Yl + ^TT-' 



m'l n' + m'l 



— \i^is{i]l + mli,m!.2]Si,\2)\ . 



In anticipation of the discussion of renormalization in Sec. bare masses mo and mi and 
bare couplings 60 and 61 have been introduced. These equations are solved numerically, with 
the first step being the use of ( |2.20| ) to eliminate 7l)2 from ( |2.19| ). Once a solution is obtained 
for one value of total spin s, the solution for the opposite spin can be computed directly 
from 



V'u(si = -1/2, Ai = Tl) = -V^i%(+1/2,±1) 



and 



V^u(si = +1/2, Ai = Tl) = +^i*t(-1/2, ±1) 
These follow from the sjTiimetries of the integral equations. 



(2.21) 



(2.22) 
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Other symmetries lead to relationships between amplitude components, which can be 
summarized as follows: 



ipi^ix, k 



_;si = 1/2, Ai = ±1) = Kfr\A\K\, \ky\) +ikyfixi{\ka:\, \ky\) 



(2.23) 



and 



i^i^ix, k±; si = -1/2, Ai = ±1) = grx^i\k^\, \ky\) + ik^kygix^{\k^\, \ky\) 



(2.24) 



where the functions fr\, fix, grx, and gix are real. The problem can then be reduced to a 
smaller matrix problem for these real functions. For Si = 1/2 we store kxfr±{kx > 0,0), 
kyfi±{0, ky > 0), kxfr±{kx > 0,ky > 0), and kyfi±{kx > 0,ky > 0). For si = —1/2 we store 
5'r±(0,0), gr±{k:i: > 0,0), ^(^±(0, /c^; > 0), gr±{kx > 0,ky > 0), and kxkygi±{kx > 0,ky > 0). 
The use of symmetry reduces the matrix storage requirement by a factor of 8. The Hermitian 
matrix of the original eigenvalue equation (|2.19| ) can be expressed as a real symmetric matrix 
in the reduced equation by using a two-component representation of complex arithmetic: 



(c + id){a + i[3) 




(2.25) 



and 




(2.26) 



The leading perturbative result |20 



is recovered by keeping only ipQs terms on the right- 
hand side of ( p.l9|) . This equation can then be immediately solved for ipis, which can 
be used to form a discrete approximation to ( |2.12| ). The approximation includes a finite 
difference approximation to the derivatives that appear in ( |2.10| ) and therefore is not simply 



a trapezoidal approximation to ( ^.12| ). 



C. Discretization errors 

Results obtained with ordinary DLCQ show an irregular dependence on the numerical 
parameters K and N±, which interferes with extrapolation to infinite resolution. The causes 
of the irregularities have been determined to be the numerical approximation of the derivative 
in the formula ( p.lO[ ) for ag and boundary effects in the numerical integrations. The error 



in the derivative can be controlled by choosing iV^ > 7 and K > 21. The bound on K is 
consistent with the resolution needed to resolve the one-photon peak in the integrand 
when = mg/lO, which is the photon mass we use. Smaller values of shift the peak 
to smaller photon momenta and would increase the lower bound on K. The shape of the 
integrand for various values of is illustrated in Fig. ^. The mass sensitivity of the 
numerical convergence rate is shown in Table |. 

The integration boundary effects are more difficult to control. These effects arise from 
use of the DLCQ grid which is incommensurate with the integration domain. At the bound- 
aries, the trapezoidal rule misses contributions beyond the last grid point; this error is not a 
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FIG. 4. The integrand for the one-photon perturbative contribution. The integrals over the 
transverse momentum have been performed and only the integral over longitudinal momentum x 
remains, as given in Eq. (2.13). 



14 



TABLE I. The DLCQ approximation to the perturbative one-photon integral in Eq. ( ^.131) . 
The values at different DLCQ resolutions K show that the convergence rate depends on the photon 
mass rriry. 



K 





0.001 


0.01 


0.1 


0.3162 


11 


0.8182 


0.8182 


0.8173 


0.7455 


0.5041 


21 


0.9048 


0.9047 


0.9025 


0.7708 


0.5079 


41 


0.9512 


0.9512 


0.9461 


0.7728 


0.5090 


81 


0.9753 


0.9752 


0.9648 


0.7731 


0.5093 


161 


0.9876 


0.9873 


0.9699 


0.7732 


0.5094 


321 


0.9938 


0.9933 


0.9703 


0.7732 


0.5094 


641 


0.9969 


0.9959 


0.9703 


0.7732 


0.5094 


1281 


0.9984 


0.9968 


0.9703 


0.7732 


0.5094 


2561 


0.9992 


0.9969 


0.9703 


0.7732 


0.5094 


CXD 


1.0000 


0.9969 


0.9703 


0.7732 


0.5094 



smooth function of the grid spacing. To overcome this error, one can replace the trapezoidal 
rule by open-closed Newton-Cotes formulas tailored specifically to the boundary [l^. Grid 
points near the boundary are then associated with unequal integration weights. The un- 
equal weights must be taken into account in normalization sums and symmetrization of the 
Hamiltonian matrix, but this is easily done. One can even consider use of Simpson's rule, al- 
though this does not appear useful in the anomalous moment calculation. The improvement 
brought by these weighting methods can be dramatic, as shown in Ref. |T^. 



III. RENORMALIZATION 

A. Mass renormalization 

Here we are interested in ultraviolet divergences associated with large Electron 
self-energy contributions, which are divergent, shift the mass and, through wave function 
renormalization, change the coupling. These induce and log A dependencies in the eigen- 
values. In the discrete truncated problem these effects depend on the Fock sector considered 
||1CI|| . For example, an electron in a Fock state for which a transition to a state with more 
photons is not allowed, perhaps due to truncation in photon number, will not experience 
any self-energy corrections. If one additional photon is allowed, but not instantaneous in- 
teractions, only single loops can occur. If two or more additional photons can appear, then 
an infinite number of overlapping loops can contribute to self-energy corrections, a truly 
nonperturbative situation. In each case, the leading divergence is removed by introduc- 
tion of counterterms associated with bare masses that are sector dependent and momentum 
dependent [|67| . 

To be specific, consider the case where there are at most two photons and only one 
electron. The Fock-state expansion can be written schematically as 
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= Vo|e) + ViIct) + ^^21677) • (3-1) 

Here i/ji and i/'2 are column vectors that contain the amphtudes for individual Fock states 
with one and two photons, respectively. The eigenvalue problem ( |2.4| ) becomes a coupled 
set of three integral equations ( p.l8|) , ( p.l9|) , and ( p.20| ), which we write more compactly as 



bi^Po + An'01 + Ai2'02 = MVi , (3.2) 
62^0 + Al^ipi + A22'02 = MV2 , 

where mo is the bare electron mass and the vectors 6| and the tensors Aij are integral 
operators obtained from Hi^c- We now require that mo be such that = m^ is an 
eigenvalue. The second and third equations can be solved for if^i/ipo and '02/V'o- Then the 
first equation yields m^. Normalization of fixes the value of ipo- 

Suppose now that this two-photon problem is embedded in some larger problem where 
one needs to know the bare mass of an electron in a Fock state that can couple to Fock 
states obtained by adding at most two photons. The same set of equations can be applied, 
with all constituents in the lowest Fock state, other than the electron, acting as spectators. 
One need only replace m^ in (|3.2|) by (mg + kjj/x and = by = (m^ + kjj/x, 
with X and k± the longitudinal momentum fraction and transverse momentum of the initial 
electron. Notice that mo is now a function of x and k±. 

This can be generalized to cases with more photons, and reduced to the case of only one 
contributing photon. Thus one obtains a mechanism for a sector-dependent, momentum- 
dependent mass renormalization that is used from the top n-photon sector down to the bare 
electron state |e). The last step automatically includes the solution of the full eigenvalue 
problem for the dressed electron state. 

For the one-photon case embedded in the two-photon problem we have 

"''"^^' ^i,(x, fcx) + bt(x, k^) ■ V2 = '"'^^' ^i.(x, k^) , (3.3) 

X X 

2 / 2 

b^uix, k^) + AV2 = "^'^ ^ i^2s . (3.4) 

X 

The second photon is a spectator. The coupling to i/'2 then induces the one-loop self-energy 
correction with this spectator present. The explicit form is obtained from (|2.19| ) and (p.20|) 



with = {ml + k"^) / X + {m'i + kjj / {1 ~ x) and with any interaction involving the spectator 



dropped. Eq. |2.20| can then be solved for -02 and the result substituted into the modified 
( ^.lij| ) to obtain. 

— V^is 21,I2i;si, Ai = 77^; ^is{n,m-Si,Xi) 3.5 

n' / K n' / K 

2, i-.K'^ml^ 5n+rn,n'/rn{l/n' -l/nf 'ipis{'nl,ml-,si,\i) 



+e?(n')frJE 



5nA 


-m,n'M {( 


m±/m - 


- n^/nf' 


+ {m^/m — n'j_/n'] 


^} -ipuin/jml; si, Ai) 


K 
n' 


[ml + (n' 




- ^ [ml 


+ [n^n/L^y] - f 


m2 + {m±7i/L±y 
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The Kronecker deltas from helicity conservation have been used to simplify the result, and 
only terms in which the second photon is a spectator have been kept. Rearrangement of the 
coefficient of ipu yields 

mlin') = ml-^el{nf)^ (3.6) 

Sn+m,n' {nieL^/Tif {l/n' - l/uf + {m±/m - n^/nf + {m±/m - n'^/n'f 
^ f K + (nl7r/L^)2] - f K + {n^-K/L^f] - f [m^ + (m^7r/L^) 

as the one-loop mass. 

If electron-positron pairs are included, the photon mass is renormalized and must be 
treated in the analogous fashion. In general, the two mass renormalizations are coupled, 
and must be carried out simultaneously. 

All of the steps in mass renormalization depend on knowing all couplings. This informa- 
tion is actually not immediately available because the couplings are to be renormalized. 

B. Coupling renormalization 

The bare coupling for the electron-photon three-point vertex depends on the initial and 



final momenta of the electron and on the sectors between which the coupling acts |T0| . The 
momentum dependence is present because the amount of momentum available constrains 
the extent to which higher order corrections can contribute. Similarly, the sector depen- 
dence makes itself felt when the number of additional particles in higher-order corrections 
is restricted. 

We fix these bare coupling functions by matching photon absorption amplitudes to the 
fundamental three-point vertex. The amplitudes are computed from the numerical eigen- 
function of the light-cone Hamiltonian. Therefore, the coupling renormalization conditions 
and the mass eigenvalue problem form a coupled set of equations that are solved iteratively. 



1. Renormalization conditions 

When vacuum polarization is absent, the bare coupling cq is related to the physical 
coupling cr by 

eo{k,kf) = ^ ' (3-7) 

where = {kf,k±i) is the initial electron momentum and kj the final momentum. The 
renormalization functions Zi{k) and Z2{k) are generalizations of the usual constants [^]. 

The wave function renormalization function Z2 is easily computed since it is the proba- 
bility of the bare electron Fock state in the dressed electron state. In the earlier notation of 
Eq. (pj.lD, we have 

Z2{k) = \tPo\\ (3.8) 
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where k is the hght-cone momentum of the dressed electron. The amphtude i/jq must be 
computed in a basis where only allowed particles appear. For example, if the vertex is the 
photon absorption process, must be computed with one less photon in the basis than in 
the basis used for Z2/. From this example one can see that the Tamm-Dancoff approximation 
has destroyed the usual Ward identity. 

The function Zi can be fixed by considering the transition amplitude Tj, for photon 
absorption by an electron at zero photon momentum. The proper part of this amplitude, 
meaning that without self-energy corrections to the legs, is required to be proportional to 
the elementary three-point no-fiip vertex Vfi when q = kj — ^ 0: 

In the limit, only kj = k^ dependence can remain.[] Numerically the limit can be taken 
by using a photon with momentum {2P^ / K^q^Ti / L_[_)] in the DLCQ limit of — * 00 and 
L_|_ — i> 00, this momentum becomes zero. Of course, if this particular state does not satisfy 
the cutoff, a state with slightly larger longitudinal momentum must be used instead. 

The full transition amplitude can be computed from solutions to the eigenvalue problem 
(^.41). Let Hq be the free light-cone Hamiltonian with physical masses. The eigenstates of 



Hq are then the asymptotic states of the electron and photon. The transition is driven by 
the interaction V = i^LC ~ Hq. Define resolvents for the free and full Hamiltonians as 

G+ = ^— - and = . ^ , (3.10) 

s + ie- Hq s + ie - Hi^c 

with s the square of the center-of-mass energy. The T matrix can be formally expressed in 
terms of these as 

G+TG+ = G+VG+ + G+VG+TG+ = g+VG+ . (3.11) 
When sandwiched between the initial and final states, this yields 

-^^Tf^^ = E ^^(^nlV^K)^ , (3.12) 
s - mj ' s - Si ^ s - s- Si 

where the |\E'„) are eigenstates of i^LC with eigenvalues and bare-electron amplitudes 
ipno- In the limit0that s becomes ml, we obtain 

Tf, = MW\^), (3.13) 



in which |\E') is the dressed electron state and i/jq = JZ2f{kf). 



^There are, of course, finite corrections that are not properly represented here. These are discussed 



in Sec. [TV . 

^This limit neglects the small photon mass. 
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FIG. 5. Representative diagram for the extraction of the proper vertex ampUtude J'P™?'^'" from 
the full amplitude Tjj. The vertical lines separate the three regions of the diagram. To the left the 
initial photon is a spectator to the dressing of the initial electron. To the right only corrections to 
the final electron line are present. The proper vertex is in the center. 



The connection between Tji and Tj^°^'^'^ is made by considering the matrix element of 
G+TG+ = G+VE^=o (G+V^)"G+. We have 

(/|G+TG+|i) = J2{f\G+V ■ ■ ■ VG+\f){f\VG+V ■ ■ ■ G+\t){t\G+V ■ ■ ■ VG+\i) . (3.14) 

The factors on the right are illustrated in Fig. ^. The second factor contains no intermediate 
I/) states and the initial photon is absorbed before \i) appears as an intermediate state. In 
the third factor, the initial photon remains a spectator throughout. The sum runs over all 
possible combinations of these forms and yields 

(/IG+TG+K) = {f\g+\f)Tfr'''{e.\g:^\e.), (3.15) 

where is the propagator for the electron in the presence of the initial photon as a spec- 
tator. In the limit where s and Sj approach ml, we obtain 



rpproper 



|2 



Tfi^— = JZI^lL^Tfr"'^^^ . (3.16) 



9 9 fi 

s — m^ s — Si s — m^ ■' s — s 
This then reduces to an expression for Tj^°^^^ 

Tfr"' = y^Tf, . (3.17) 

Thus the solution of the eigenvalue problem for only one state can be used to compute Zi. 
Full diagonalization of i^LC is not needed. 
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2. Application of renormalization conditions 



Because Zi is needed in the construction of i/LC? the eigenvalue problem and the renor- 
malization conditions must be solved simultaneously. This leads to an iterative procedure 
that begins with an initial guess for the bare coupling functions. One then computes bare 
masses and new bare couplings. The process is repeated until convergence is attained. This 
must be done from the top sector down; the bare masses in any one sector and the bare 
couplings between any two depend only on the sectors above, the ones with more photons. 
The structure of the Hamiltonian matrix can then be determined once and for all at these 
levels and then used in the determination of the structure at the levels further below. 

When the Fock basis is limited to no more than one photon, and instantaneous interac- 
tions are neglected, the renormalization conditions are quite simple. We have from ( |3.13| ) 



and from this, with ( |3.9D and (|3.17| ), 



Z\ — Z2fZ2i 



Z2fVfi 



'2i 



1 



(3.18) 



(3.19) 



where the last equality follows from the unavailability of any state that can correct the initial 
electron line when a photon spectator is present. The bare charge is then given by 



ei 



ZiCr 



Z2fZ2i ^0 



(3.20) 



The subscript of 1 corresponds to use in couplings between one and two-photon states. 

We now consider the solution of the problem in the case of a basis with no more than two 
photons. Eq. (|3.20|) provides the solution for the bare coupling between one and two-photon 
states and, through spectator dependence of '?/'o, makes ei a function of the final electron 
momentum. We then need to consider the bare coupling between the bare electron and the 
one-photon states. On substitution of Eqs. ( |3.9| ) and ( |3.17| ), Eq. ( ^.71 ) becomes 



eo 



Z2fZ2i7T^^R ■ 



(3.21) 



fi 



This is a nonlinear equation for eo because Tfi has a complicated dependence on this bare 
charge. To make this dependence explicit, we first use the fact that 1^^) is an eigenstate of 
Hq + V to reduce ( |3.13| ) to the form 



T 



fi 



ml + TcyL] 
1-2/K 



m2 + nyL] 



2/K 



iJori{l-2/K,xn/L^) 



(3.22) 



The amplitude ipi satisfies the middle equation of (pl2|), which can be written as 

biV'o + [A' + elcc^] ■ -01 = m^Vi , (3.23) 
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where A' is an effective interaction obtained by integrating out the amphtude and the 
CqCC^ term is the finite instantaneous fermion interaction. The couphng to ipo contains eg 
as a simple factor. We extract this in the following definitions: 



b[ = bi/eo and ip[ = -01/(60^0) • 



(3.24) 



The scaled amplitude can then be obtained as a formal solution to Eq. ( p.23|) that shows 
all Co dependence^ 



B 



D 



with 



B 



{A' 



m 



h\ and D = {A - mlY^c 



The amplitude for two-photon states is then given by 

where i? is a rectangular matrix independent of cq. 
The discrete normalization condition is 



(3.25) 



(3.26) 



(3.27) 



1 = l^oP + + |^2p = |^oP[l + elm' + eglV'^ 



This yields 



which can be used with ( p.22|) and (|3.21| ) to obtain 



(3.28) 



(3.29) 



Co 



ip'{{l-2/K, X7i/Lj 



1-2/K 



2/K 



(3.30) 



where V = V/eQ is independent of Cq. The phases of ip[ and Vj^ are such that the right- 
hand side is real, as it must be. The remaining implicit dependence on cq is in i/?'^, which 
is given by ( |3.25| ), and in t/>2 = Rtj^'i, with R independent of cq. Notice that B and D are 
independent of cq and need to be computed only once. The equation for eo is best solved 
iteratively after it is squared to eliminate the square root on the right hand side. 

A real solution exists only for a finite range of the physical coupling cr. This is an 
artifact of the Tamm-Dancoff truncation and the consequent failure of the Ward identity. 
The value of the critical coupling e^^'*, the upper limit of the allowed range, can be found 
by studying the eo — > 00 limit. In this limit we find 



^Notice that A' is independent of eg, as can be seen from (2.1£) and ( p.20| ). 
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a 



crit 
R 



47r 47rZ2.r/2(|^;|2+|V>'2|2 



_ m2 + ttVL^ 

1 1 



m: 



2/K 



2/K 



(3.31) 



with calculable as S — (c''" ■ B / ■ D)D. Values of the critical coupling are plotted in 
Fig. 1^. The change from a basis with no more than one photon to a basis with no more than 
two is quite small. To stay within the limit imposed by this result we will use a physical 
coupling of aR = 0.1. 



C. Infrared singularities 

A nonzero photon mass is used to eliminate the usual infrared singularities. Because 
the calculation deals with a charged system, there would otherwise be considerable difficulty 
with soft photons ||68|. On the light cone, there are other singularities not removed by 
the photon mass. They are associated with contributions that involve zero longitudinal 
momentum. 

The fundamental four-point vertices can be infrared singular, in the limit of zero longi- 
tudinal momentum for the instantaneous fermion. They must be allowed to cancel against 
iterations of the three-point vertices which are also singular. This constrains the bare cou- 
plings in the four-point vertices to forms derived from the three-point coupling cq. The 
pairs of diagrams are shown in Fig. |^. The first pair does not actually involve a singularity; 
however, we do match the four-point coupling to the iterated three-point coupling. The 
second requires basis states with two photons, which are available in the calculation. The 
third pair requires the presence of electron-positron pairs in the basis, or an effective inter- 
action in the Hamiltonian. Neither is included at present, and therefore this piece of the 
instantaneous interaction must also be excluded from the Hamiltonian. Detection of < 
in the instantaneous interaction can be easily done to exclude this graph. 

Other infrared singularities are associated with the emission and absorption of real pho- 



tons with longitudinal momentum near zero In perturbative calculations these are 

regulated by the Mandelstam-Leibbrandt prescription ||70[. Viewed in a;"'"-ordered pertur- 
bation theory, each intermediate state contributes a denominator in which the light-cone 
energy k~ = (m^ + k'^)/k^ of the photon becomes large, and each vertex can contain a fac- 
tor of (k~^)~'^^'^. If k~^ is separately regulated, with some lower cutoff e, graphs with multiple 
photons will contribute powers of loge or even e~^. For a Tamm-Dancoff approximation 
to a charged system, these cannot be expected to cancel. The choice of the invariant mass 
cutoff ( p.ll|) instead couples the regulation of small k~^ to that of large k± . The combination 
prevents the small k~^ region of integration from making large contributions except in cases 
where there are already ultraviolet transverse divergences. These spurious infrared infinities 
are then handled by the mass and coupling renormalization discussed in this section. 



D. Four-point graphs 

There remains a logarithmic divergence associated with four-point graphs of the sort 
illustrated in Fig. ^(a) and (b). If all graphs of this order are included in a perturbative 
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FIG. 6. The critical coupling ag'* as a function of cutoff A. This is the upper limit of 
the range allowed for the physical coupling; beyond this value there is no solution for the bare 

coupling, as a consequence of the truncations of the theory. The one-photon values are obtained 
from extrapolations using a^'* + ai/K + a2/K'^ + i>i/N±_ + 62/-^! + Cii/{KN_[_); the error bars 
represent the difference between this fit and one without the bi/N± term. The two-photon values 
come from a fit to ag'* + ai/K + a2/K^ + bi/N± + b2/N]_; the error bars are obtained from a fit 
without the quadratic terms. For all the photon mass is rrij = me/10. 
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FIG. 7. Instantaneous fermion interactions paired with corresponding iterated three-point 
interactions. The vertical dashed lines indicate the intermediate states of the iterated interactions. 
In (b) the longitudinal momentum of the instantaneous fermion is positive, and in (c) it is negative; 
this represents a separation of the crossed-photon graph into two pieces. 

calculation, the logarithms cancel. However, the Tamm-Dancoff truncation of the present 
calculation excludes some graphs, such as the one shown in Fig. |](c), and the cancellation 
can no longer take place. In a nonperturbative calculation one must include the equivalent 
of diagrams with an arbitrary number of interlocked loops, such as Fig. ^(b), which are also 
logarithmically divergent. The needed counterterm is of the form X{pt,p~f) log A but cannot 
be found analytically without summing all the infinite chains of interlocked loops. 

One way to approach the construction of X{pt,p~f) is to fit Compton amplitudes to 
data [0] . This will require development of techniques to describe scattering processes within 
DLCQ. A generalization of earlier work ||45|j7l[| on the inversion of the full Greens function 
may be useful as a means for computing the T matrix and thus scattering amplitudes. We 
do not consider this further here. The results presented in Sec. do still contain the log A 
divergence. 



IV. FINITE CORRECTIONS 

A. Photon zero modes 

As applied to QED, DLCQ requires the use of periodic boundary conditions for the 
photon field. This is because photons couple to fermion bilinears, which are automatically 
periodic, even if the preferred antiperiodic boundary conditions are used for the fermions. 
For fields periodic in the longitudinal direction x~ , there are contributions from zero modes 
1 26| , |57| -p0| , modes independent of x~ that correspond to zero longitudinal momentum. As 



shown by Pauli and Kalloniatis [^], these modes prevent the choice of ordinary light-cone 



gauge because the zero-mode piece of cannot be gauged away. Instead this piece must 
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FIG. 8. Logarithmically divergent four-point graphs. Those with no more than two photons in 
any intermediate state, (a) and (b), are included in the present calculation. The nonperturbative 
nature of the calculation implies that (b) can have an arbitrary number of loops. Diagram (c), 
which contains three photons in flight is not included. 
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satisfy a constraint equation. In fact, careful application of DLCQ to most bosonic theories 
will result in constraint equations that relate the zero-mode contribution to the normal- 
mode operators in a nonlinear, nontrivial way. For QED there are zero modes in all the 
components of the photon field, for which constraint equations must be solved. What is 
more, the constraint equation for the dependent piece of the fermion field, which is easily 
solved in light-cone gauge in the continuum, becomes coupled to the zero-mode constraint 
equations. A formulation of the coupled system of constraints has been given by Kalloniatis 



and Robertson |^6[. Extension of these constraints to include a nonzero photon mass is 
st raight forward . 

The constraint equations are difficult to solve, even in simpler theories ||6^. This is 
partly because they couple states with different P+ and require study of convergence as a 
function of some cutoff. The difficulty is also due to the need for an ultraviolet cutoff 
and renormalization of masses and couplings. Because the renormalization is formulated 
in terms of solutions to the mass eigenvalue problem and because the Hamiltonian cannot 
be formed until the zero-mode contribution is known, the problem expands to a very large 
nonlinear system of simultaneous equations. As a result of these difficulties, the calculations 
discussed here do not include zero modes. However, some progress has been made recently 
by Kalloniatis in the solution of constraint equations for SU(2) Yang-Mills in 1 + 1 
dimensions coupled to massive adjoint scalars. 

For theories such as QED where symmetry breaking effects are not expected, solution 
of constraint equations may not be necessary. One can instead treat the end-point behav- 
ior of photon amplitudes in a manner similar to that of the "ladder relations" studied by 
Antonuccio and Dalley [^. Behavior of amplitudes at small longitudinal momentum, as 
extracted from the integral equations, can be used to construct effective interactions that 
include zero modes to leading order in 1/K. This is equivalent to the approach used in PB 
where the behavior of the exchange kernel was studied in a scalar theory to determine the 
effective interaction |73|. To keep zero-mode terms to higher order in 1/K would actually 



be inconsistent with DLCQ's neglect of higher order non-zero-mode terms. In the work of 
Ref. [35 1 inclusion of the zero-mode contributions 0{1/K) did improve convergence. 

The whole issue may actually be moot when the invariant mass cutoff (|2.11|) is used. This 
cutoff explicitly excludes contributions from states with zero longitudinal momentum. The 
meaning of this exclusion for nondynamical fields is unclear. The calculations that showed 
zero modes to be useful for convergence ||3^ did not employ the invariant mass cutoff. New 
calculations need to be carried out specifically to study the effect of cutoff choice on the 
importance of zero modes. 



B. Restoration of symmetries 

The use of light-cone coordinates, combined with the Tamm-Dancoff truncation in parti- 



cle number and the invariant mass cutoff, explicitly break symmetries of the theory ||10[ . In 



particular, rotational symmetry about the transverse axes is broken because the associated 
operators involve the interaction and therefore change particle number. The change in par- 
ticle number cannot be accommodated in field theory without allowing an infinite number 
of particles. 
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Restoration of such symmetries can be accomplished by the addition of finite counter- 
terms to the Hamiltonian 0,^] including adjustment of the "vertex mass," which appears 
in the spin- flip vertex, relative to the "kinetic mass" |]73[. The ambiguities associated with 
the infinite counterterms allow such finite terms to exist [|l^. Restoration of symmetries 



is then viewed as a source of conditions by which these finite parts can be determined. 
In practice, this might involve study of processes such as Compton scattering or 
electron-electron scattering. 

Given the Tamm-Dancoff truncation, an alternative is to view the eigenvalue problem as 
a few-body problem [[7^] for which the correct effective Hamiltonian and the generators of 
translations, rotations, and boosts must satisfy the usual Poincare algebra [|7|]. The effective 
operators might be constructed by adding minimal finite corrections to their field-theoretic 
forms. The finite corrections are determined by the requirement that the Poincare algebra 
be satisfied. 

For the results presented here, no attempt has been made to include these finite correc- 
tions. 



V. RESULTS 

An accurate DLCQ calculation for a basis with at most one photon can be easily done 
when instantaneous interactions are neglected. The accuracy can be verified directly because 
the integrals that yield Oe can be performed analytically In the limit of infinite cutoff 



and zero photon mass this reproduces the Schwinger UTSf result of aij/27r. The only coupling 



renormalization is a trivial wave function renormalization. The DLCQ result at various cutoff 
values is shown in Fig. ^ for a photon mass of mg/lO. Weighting methods are a critical 
part of the calculation. 

Calculations with a basis that includes at most two photons have been done at five 
different values of the cutoff for a couphng of = 0.1. They are also shown in Fig. ^. The 
two-photon contribution adds approximately 40%. This is much larger than the order of 
magnitude (a/vr or 3%) that one would expect. It is also opposite in sign to the Sommerfield- 
Petermann contribution of — 0.328(Q;/7r)^ to the anomalous moment. We attribute this 
large difference to the absence of Z graphs. 

The basis sizes involved are on the order of 1 to 4 million, which translates to solution of 
linear systems with 4000 to 10,000 variables once the two-photon states are integrated out 
and symmetries of the one-photon states are used. 

The values obtained from DLCQ were extrapolated to i^' = oo and Nj_ = oo hj fits to 
Ge + a/K + b/N± + c/ {KN±). Exclusion of the last term provided an estimate of the error 
in the fit, which is reflected in the error bars in Fig. |^. The values of K ranged from 21 to 
31 and those of N± from 7 to 9, 10 or 11, depending on matrix size limitations. 

The time required for an extrapolated value at a fixed cutoff is roughly 10 hours on a Cray 
X-MP, using less than 32 million words of memory. This seems quite competitive with older 
lattice methods, where a quenched QCD calculation of heavy-light meson wave functions 
80| required 300 hours on a CM-200 [^, but does not yet match the effort required with 



the latest methods for which calculation of the B meson magnetic form factor might 



require 50 hours on a good PC |82 
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FIG. 9. The ratio of the renormahzed anomalous moment to the exact one-photon pertur- 
bative result af at the same photon mass m^, as a function of the invariant-mass cutoff A. Here 
= m-e/lO, and an = 0.1. The Fock space is truncated to include at most one or two photons. 



The DLCQ results are extrapolated from calculations done with K 
for one photon and A^^ = 7 to 9, 10, or 11 for two photons. 



21 to 31, and A^_l = 8 to 15 



VI. SUMMARY 



The nonperturbative calculation of the anomalous moment of the electron Og, besides be- 
ing of intrinsic interest itself, exposes many important issues for nonperturbative calculations 
within gauge theories which occur in the context of a truncated Fock space. These include 
nonperturbative mass and coupling renormalization, control of spurious infrared singulari- 
ties, determination of zero-mode contributions, and the construction of finite counterterms 
which restore symmetries. Each of these has been addressed in the preceding sections, and 
the first two have been incorporated into a DLCQ calculation where as many as two photons 
are included in the basis. 

We have presented results for Og computed in a light-cone-gauge Fock space truncated 
to include one bare electron and at most two photons; i.e., up to two photons in flight. 
The calculational scheme uses an invariant mass cutoff, discretized light-cone quantization 
(DLCQ), a Tamm-Dancoff truncation of the Fock space, and a photon mass regulator. We 
have utilized new weighting methods which greatly improve convergence to the continuum 
within DLCQ. A large renormalized coupling strength an = 0.1 is then used to make the 
nonperturbative effects in the electron anomalous moment from the one-electron, two-photon 
Fock state sector numerically detectable. Results are given in Fig. ^. 

The disagreement between these results and what one would expect from perturbation 
theory at order indicates that the effect of Z graphs needs to be included in a systematic 
way. This can be done as an effective interaction, to avoid expansion of the Fock basis to 



28 



include pair states. The corresponding piece of the instantaneous fermion interaction, as 
depicted in Fig. |^(c), must then also be included to maintain an infrared cancellation. 
Further progress in computing the electron moment will require: 

1. New counterterms: One piece of the infinite renormalization is missing in the cal- 
culation. As discussed in Section [111 D|, it requires a new nonperturbative counterterm for 



the logarithmic divergences present in diagrams of the type shown in Fig. |g(a) and (b). 
The divergence arises because the Tamm-Dancoff truncation prevents certain cancellations. 
Construction of the counterterm will likely require analysis of scattering processes. 

2. Zero modes in DLCQ: Before full consideration of photon zero modes is undertaken, 
we recommend renewed study of zero modes in a scalar theory where the constraint equation 



can be solved exactly [^. This may show that, when the invariant mass cutoff (|2.11| ) is used, 
zero modes do not make a significant numerical contribution. If instead there is an important 
contribution, it should be computed only to leading order in the numerical resolution, to be 
consistent with the level of approximation used in the basic DLCQ approach. 

3. Use of symmetries in DLCQ renormalization: The restoration of symmetries should 
then complete construction of the light-cone Hamiltonian. One can normalize to specific 
physical processes [ff6H7| or take an abstract approach based on the algebra of the Poincare 
generators [ [77| , [75 |. 



4. Higher Fock States: Once the two-photon calculation is fully under control, the 
addition of eee^ states can be considered. This will require analysis of photon mass and 
wave function renormalization. 

Many of the complications of the light-cone Fock state analysis presented here can be 
traced to the complexity of sector-dependent renormalization. Given such complications 



the newly developed alternative of Pauli-Villars regularization [|T^ may be the preferred 
approach. Within such a scheme, it is also likely that the limitation to a small number of 
photons can be relaxed. 

The analysis presented here is the first step in a systematic program to compute physical 
quantities in gauge theory systematically utilizing a light-cone Fock expansion. It will also 
be interesting to use these methods and the present knowledge of the dressed-electron state 
in QED in order to systematically construct the neutral positronium state as a composite 
of a dressed electron and positron. Such an analysis can serve as the prototype for systemic 
nonperturbative construction of colorless bound state hadrons in QCD. 
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